Chapter One 


Introduction 


This book represents a modern treatment of classical control theory of continuous- 
and discrete-time linear systems. Classical control theory originated in the fifties 
and attained maturity in the sixties and seventies. During that time control theory 
and its applications were among the most challenging and interesting scientific 
and engineering areas. The success of the space program and the aircraft industry 
was heavily based on the power of classical control theory. 

The rapid scientific development between 1960 and 1990 brought a tremen¬ 
dous number of new scientific results. Just within electrical engineering, we 
have witnessed the real explosion of the computer industry in the middle of the 
eighties, and the rapid development of signal processing, parallel computing, 
neural networks, and wireless communication theory and practice at the begin¬ 
ning of the nineties. In the years to come many scientific areas will evolve 
around vastly enhanced computers with the ability to solve by virtually brute 
force very complex problems, and many new scientific areas will open in that 
direction. The already established “information superhighway” is maybe just 
a synonym for the numerous possibilities for “informational breakthrough” in 
almost all scientific and engineering areas with the use of modern computers. 
Neural networks—dynamic systems able to process information through large 
number of inputs and outputs—will become specialized “dynamic” computers 
for solving specialized problems. 

Where is the place of classical (and modern) control theory in contemporary 
scientific, industrial, and educational life? First of all, classical control theory 
values have to be preserved, properly placed, and incorporated into modern scien¬ 
tific knowledge of the nineties. Control theory will not get as much attention and 
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recognition as it used to enjoy in the past. However, control theory is concerned 
with dynamic systems, and dynamics is present, and will be increasingly present, 
in almost all scientific and engineering disciplines. Even computers connected 
into networks can be studied as dynamic systems. Communication networks have 
long been recognized as dynamic systems, but their models are too complex to 
be studied without the use of powerful computers. Traffic highways of the fu¬ 
ture are already the subject of broad scale research as dynamic systems and an 
intensive search for the best optimal control of networks of highways is under¬ 
way. Robotics, aerospace, chemical, and automotive industries are producing 
every day new and challenging models of dynamic systems which have to be 
optimized and controlled. Thus, there is plenty of room for further development 
of control theory applications, both behind or together with the “informational 
power” of modern computers. 

Control theory must preserve its old values and incorporate them into modern 
scientific trends, which will be based on the already developed fast and reliable 
packages for scientific numerical computations, symbolic computations, and 
computer graphics. One of them, MATLAB, is already gaining broad recognition 
from the scientific community and academia. It represents an expert system 
for many control/system oriented problems and it is widely used in industry 
and academia either to solve new problems or to demonstrate the present state 
of scientific knowledge in control theory and its applications. The MATLAB 
package will be extensively used throughout of this book to solve many control 
theory problems and allow deeper understanding and analysis of problems that 
would not otherwise be solvable using only pen and paper. 

Most contemporary control textbooks originated in the sixties or have kept 
the structure of the textbooks written in the sixties with a lot of emphasis 
on frequency domain techniques and a strong distinction between continuous- 
and discrete-time domains. At the present time, all undergraduate students in 
electrical engineering are exposed to discrete-time systems in their junior year 
while studying linear systems and signals and digital signal processing courses 
so that parallel treatment of continuous- and discrete-time systems saves time 
and space. The time domain techniques for system/control analysis and design 
are computationally more powerful than the frequency domain techniques. The 
time domain techniques are heavily based on differential/difference equations and 
lineal - algebra, which are very well developed areas of applied mathematics, for 
which efficient numerical methods and computer packages exist. In addition, the 
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state space time domain method, to be presented in Chapter 3, is much more 
convenient for describing and studying higher-order systems than the frequency 
domain method. Modern scientific problems to be addressed in the future will 
very often be of high dimensions. 

In this book, the reader will find parallel treatment of continuous- and 
discrete-time systems with emphasis on continuous-time control systems and on 
time domain techniques (state space method) for analysis and design of linear 
control systems. However, all fundamental concepts known from the frequency 
domain approach will be presented in the book. Our goal is to present the 
essence, the fundamental concepts, of classic control theory—something that 
will be valuable and applicable for modern dynamic control systems. 

The reader will find that some control concepts and techniques for discrete¬ 
time control systems are not fully explained in this book. The main reason 
for this omission is that those “untreated topics” can be simply obtained by 
extending the presented concepts and techniques given in detail for continuous¬ 
time control systems. Readers particularly interested in discrete-time control 
systems are referred to the specialized books on that topic (e.g. Ogata, 1987; 
Franklin et ai, 1990; Kuo, 1992; Phillips and Nagle, 1995). Instructors who 
are not enthusiastic about the simultaneous presentation of both continuous- and 
discrete-time control systems can completely omit the “discrete-time parts” of 
this book and give only continuous-time treatment of control systems. This book 
contains an introduction to discrete-time systems that naturally follows from their 
continuous-time counterparts, which historically are first considered, and which 
physically represent models of real-world systems. 

Having in mind that this textbook will be used at a time when control theory 
is not at its peak, and is merging with other scientific fields dealing with dynamic 
systems, we have divided this book into two independent parts. In Chapters 2-5 
we present fundamental control theory methods and concepts : transfer function 
method, state space method, system controllability and observability concepts, 
and system stability. In the next four chapters, we mostly deal with applications 
so that techniques useful for design of control systems are considered. In Chapter 
10, an overview of modern control areas is given. A description of the topics 
considered in the introductory chapter of this book is given in the next paragraph. 

Chapter Objectives 

In the first chapter of this book, we introduce continuous- and discrete-time 
invariant linear control systems, and indicate the difference between open-loop 
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and closed-loop (feedback) control. The two main techniques in control system 
analysis and design, i.e. state space and transfer function methods, arc briefly 
discussed. Modeling of dynamic systems and linearization of nonlinear control 
systems are presented in detail. A few real-world control systems are given 
in order to demonstrate the system modeling and linearization. Several other 
models of real-world dynamic control systems will be considered in the following 
chapters of this book. In the concluding sections, we outline the book’s structure 
and organization, and indicate the use of MATLAB and its CONTROL and 
SIMULINK toolboxes as teaching tools in computer control system analysis and 
design. 

1.1 Continuous and Discrete Control Systems 

Real-world systems are either static or dynamic. Static systems are represented 
by algebraic equations, and since not too many real physical systems are static 
they are of no interest to control engineers. Dynamic systems are described either 
by differential/difference equations (also known as systems with concentrated or 
lumped parameters) or by partial differential equations (known as systems with 
distributed parameters). Distributed parameter control systems are very hard 
to study from the control theory point of view since their analysis is based 
on very advanced mathematics, and hence will not be considered in this book. 
At some schools distributed parameter control systems are taught as advanced 
graduate courses. Thus, we will pay attention to concentrated parameter control 
systems, i.e. dynamic systems described by differential/difference equations. It 
is important to point out that many real physical systems belong to the category 
of concentrated parameter control systems and a large number of them will be 
encountered in this book. 

Consider, for example, dynamic systems represented by scalar differen¬ 
tial/difference equations 

x{t) = fc(x(t)), x(t 0 ) = x 0 (1.1) 

x(k + 1) = f d (x(k)), x(k 0 ) = x 0 (1.2) 

where t stands for continuous-time, k represents discrete-time, subscript c in¬ 
dicates continuous-time and subscript d is used for discrete-time functions. By 
solving these equations we learn about the system’s evolution in time (system 
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response). If the system is under the influence of some external forces, or if 
we intend to change the system response intentionally by adding some external 
forces, then the corresponding system is represented by the so-called controlled 
differential/difference equation, that is 


x(t) = fc(x(t), u(t)), x(t 0 ) = x 0 (1.3) 

x(k + 1) = fd(x(k), u(k)), x(k 0 ) = x 0 (1.4) 

where u(t) and u(k) play the role of control variables. By changing the control 
variable we hope that the system behavior can be changed in the desired direction, 
in other words, we intend to use the control variable such that the system response 
has the desired specifications. When we are able to achieve this goal, we are 
actually controlling the system behavior. 

The general control problem can be formulated as follows: find the control 
variable such that the solution of a controlled differential/difference equation has 
some prespecified characteristics. This is a quite general definition. In order to be 
more specific, we have to precisely define the class of systems for which we are 
be able to solve the general control problem. Note that the differential/difference 
equations defined in (1.1)—(1.4) are nonlinear. In general, it is hard to deal with 
nonlinear systems. Nonlinear control systems are studied at the graduate level. 
In this undergraduate control course, we will study only linear time invariant 
control systems. 

Continuous- and discrete-time linear time invariant dynamic systems are de¬ 
scribed, respectively, by linear differential and difference equations with constant 
coefficients. Mathematical models of such systems having one input and one 
output are given by 


d n y(t) 
dt n ' 


ttn—l 


d n ~ 1 y(t) 

dt 71 - 1 


+ ••• + 01 


dy(t) 

dt 


+ a 0 y(t) 


u(t) 


(1.5) 


and 


y(k + n) + (in-tyik + n - 1) +-b ayyik + 1) + a 0 y(k) = u(k) (1.6) 

where n is the order of the system, y is the system output and u is the external 
forcing function representing the system input. In addition to the external forcing 
function the system is also driven by its internal forces coming from the system 
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initial conditions. For continuous-time systems, the initial conditions are specified 
by known values of the system output derivatives up to the order of n — 1 at 
the initial time t 0 , that is 


,. , dy(t 0 ) d n 1 y(to) 

){1oh ~di~' 


(1.7) 


In the discrete-time domain the initial conditions are specified by 


y(ko), y(k 0 + 1), •••, y(k 0 + n - 1) (1.8) 


It is interesting to point out that in the discrete-time domain the initial conditions 
carry information about the evolution of the system output in time from k 0 to 
/,'d + n — 1. In this book, we study only time invariant continuous and discrete 
systems for which the coefficients a;, i = 0,1,...,?? — 1, arc constants. A block 
diagram representation of such a system is given in Figure 1.1. 


u(t) 

Input 


> 


System 


y(t) 

-1 

Output 


Figure 1.1: Input-output block diagram of a system 


In general, the input function can be differentiated by the system so that 
the more general descriptions of time invariant continuous and discrete systems 
are given by 


d n y(t) d^y^t) 

+ a n-1 - 


dt r 


dt 71 - 1 


ar-rr^ + a 0 y(t) 


dt 


d m u{t) d^u(t) 

— ”m. , . T !J m — l 


dt r 


dt 


m — 1 


+ - + b ou(t) 


(1.9) 
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and 

y(k + ??) + a n -iy(k + n - 1) +- \- aiy(k + 1) + a 0 y(k) 

( 1 . 10 ) 

= b m u(k + ???) + b m _iu(k + m — 1) + • • • + biu(k + 1) + boii(k) 

where all coefficients a;,? = 1,2,...,??, and bj,j = 0,1,...,???, are constant. 

The problem of obtaining differential (difference) equations that describe 
dynamics of real physical systems is known as mathematical modeling. In 
Sections 1.4 and 1.5 this problem will be addressed in detail and mathematical 
models for several real physical systems will be derived. 

Basic Linear Control Problem 

In summary, we outline the basic problem of the linear control theory. The 
problem of finding the system response for the given input function u(t) or 
u(k) is basically the straightforward problem of solving the corresponding linear 
differential or difference equation, (1.9) or (1.10). This problem can be solved 
by using standard knowledge from mathematical theory of linear differential 
and/or difference equations. However, the linear control problem is much more 
challenging, namely the input function u(t) has to be found such that the linear 
system response has the desired behavior. A simplified version of the above 
basic lineal - control problem will be defined in Chapter 6 using the notion of 
system feedback. The basic linear control problem can be studied either in the 
time domain (state space approach) or in the frequency domain (transfer function 
approach). These two approaches will be presented in Chapters 2 and 3. 


1.2 Open-Loop and Closed-Loop Control Systems 

It seems that if an input function can be found such that the corresponding system 
has the desired response, then the control problem of interest is solved. This is 
true, but is it all that we need? Assume that u(t) is such a function, which 
is apparently a function of time. Imagine that due to parameter variations or 
due to aging of the system components the system model is a little bit different 
than the original one or even worse that the coefficients a /. i = 0,1,2,...,?? — 1; 

hj.j = 0. 1.2. m. in equation (1.9) are not very precisely known. Then 

the function ??(/), given as a precomputed time function, might not produce 
a satisfactory solution (especially in the long run). One may try to solve the 
problem again and get a new expression for ??(/) at the expense of additional 
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computation, which is fine if the system parameters are exactly known, but 
this approach will not bring any improvement in the case when the system 
coefficients, obtained either analytically or experimentally, are known only with 
certain accuracy. This precomputed time function u(t) (or u(k) in the discrete¬ 
time domain), which solves the control problem, is known as the open-loop 
control. 

Imagine now that one is able, in an attempt to solve the basic linear control 
problem, to obtain the desired input function as a function of the system desired 
response, or even more precisely as a function of some essential system variables 
that completely determine the system dynamics. These essential system variables 
are called the state space variables. It is very natural to assume that for a system 
of order n, a collection of n such variables exist. Denote the state space variables 
by xi(t), x- 2 (t), ...,x n (t). These state variables at any given time represent the 
actual state of the system. Even if some parameters of the system are changing 
in time or even if some coefficients in (1.9) are not precisely known, the state 
variables z\(t ), x- 2 (t ), ...,x n (t ) will reflect exactly the state of the system at any 
given time. The question now is: can we get the required control variable (system 
input) as a function of the state space variables? If the answer is yes, then the 
existence of such a u(x(/)) indicates the existence of the so-called state feedback 
control. In some cases it is impossible to find the feedback control, but for 
the lineal - time invariant systems studied in this book, linear feedback control 
always exists. The linear feedback control is a linear function of the state space 
variables, that is 


«(x) = iq.x-i + F 2 x 2 H-b F n x n = Fx 

F = [F u F 2 ,.^F n ] (1.11) 

X = [.X’l, .X’2, ...,X n ] T 


where the coefficient matrix F is the feedback gain and the vector x is known 
as the state space vector. It is sometimes desirable (and possible) to achieve 
the goal by using instead of the state feedback control m xi /! i, the so-called 
output feedback control given by u(y(t)) = u(y(x(t))), which in general does 
not contain all state variables, but only a certain collection of them. In Section 
3.1 we will learn how to relate the output and state space variables. 
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Open-loop control, state feedback control, and output feedback control are 
represented schematically in Figure 1.2. 



(a) 


(b) 


(c) 


Figure 1.2: Open-loop (a), state feedback (b), and output feedback (c) controls 

The system represented in Figure 1.1 and given in formulas (1.5)—(1.10) 
has only one input u and one output y. Such systems are known as single - 
input single-output systems. In general, systems have several inputs and several 
outputs, say r inputs and p outputs. In that case, we have 

u = [u 1 ,u 2 ,...,u r ] T , y = [yi,y 2 ,...,y p ] T (1.12) 

and the matrix F is of dimension r X n. These systems are known as multi-input 
multi-output systems. They are also called multivariable control systems. A block 
diagram for a multi-input multi-output system is represented in Figure 1.3. 

Feedback control is almost always desirable as a solution to the general 
control problem, and only in rare cases and in the cases when it is impossible to 
find the feedback control has one to stick with open-loop control. Throughout of 
book we will see and discuss many advantages of feedback control. The main 
role of feedback is to stabilize the system under consideration. The feedback 
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also reduces the effect of uncertainties in the system model. In addition, it 
efficiently handles system parameter changes and external disturbances attacking 
the system, by simply reducing the system output sensitivity to all of these 
undesired phenomena. Most of these good feedback features will be analytically 
justified in the follow-up chapters. 


Ml 


_w Vi 




Mo w 


_ p.y2 

^ - w 

System 


Ilf- -^ 





-► y ‘ 


Figure 1.3: Block diagram of a multi-input multi-output system 

However, one has to be careful when applying feedback since it changes the 
structure of the original system. Consider the system described by equation (1.3) 
under the state feedback control, that is 

x(t) = f c (x(t),u(x(t))) = Fciz :(/)), x(t 0 ) = x 0 (1.13) 

Thus, after the feedback is applied, a new dynamic system is obtained, in other 
words, for different values of u(x(t)) we will get different expressions for the 
right-hand side of (1.13) so that equation (1.13) will represent different dynamic 
systems. 


1.3 State Space and Transfer Functions 

In analysis and design of linear time invariant systems two major approaches 
are available: the time domain state space approach and the frequency domain 
transfer function approach. Both approaches will be considered in detail in 
Chapters 2 and 3. Electrical engineering students are partially familial - with 
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these methods and the corresponding terminology from undergraduate courses 
on lineal - systems and digital signal processing. 

In this section only the main definitions are stated. This is necessary to 
maintain the continuum of presentation and to get mathematical models of some 
real physical systems in state space form. In addition, by being familiar with the 
notion of the state space, we will be able to present the linearization of nonlinear 
systems in the most general form. 

In the previous section it is indicated that for an nth order system there 
are n essential state variables, the so-called state space variables, which form 
the corresponding state space vector x. For linear time invariant systems, the 
vector of state space variables satisfies the linear time invariant vector differential 
(difference) equation known as the state space form. The state space equations 
of linear systems are defined by 


x(t) = Ax(t) + Bu(t), x(0) = x o 
y (t) = Cx(t) + Du(t) 


(1.14) 


x(& + 1) = Ax(t) + Bu(t), x 
y (k) = Cx(fc) + Du(A’) 

In view of the discussion in the previous section, the introduced constant matrices 
have the dimensions A nX \B nXr , C pXn , D pXr . Of course, the given matrices 
for continuous and discrete systems have different entries. In formulas (1.14) and 
(1.15) the same notation is kept for both continuous and discrete systems for the 
sake of simplicity. In the case when a discrete-time system is obtained by sam¬ 
pling a continuous-time system, we will emphasize the corresponding difference 
by adding a subscript d for the discrete-time quantities, e.g. A,/. B,/. C,/. D,/. 

The system transfer function for a single-input single-output time invariant 
continuous system is defined as the ratio of the Laplace transform of the system 
output over the Laplace transform of the system input assuming that all initial con¬ 
ditions are zero. This definition implies that the transfer function corresponding 
to (1.9) is 


( 0 ) = 


X n 


(1.15) 


b m s m + b m -\s m 1 + • • • + b\s + bo 
s n + -|-1- a x s + a 0 


(1.16) 


Polynomial exponents in transfer functions of real physical systems always satisfy 
n > in. In the discrete-time domain, the 2-transform takes the role of the Laplace 
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transform so that the discrete-time transform function is given by 

ru . b m z m + b m -\z m 1 + • • • + b\z + bo 

G(z) = -=- 

~ n + a n -iz n 1 +- \-c^z + a 0 


(1.17) 


For multi-input multi-output systems with r inputs and p outputs, transfer func¬ 
tions are matrices of order p x r whose entries are the corresponding transfer 
functions from the ith system input to the jth system output, i = l,2,...,r; 
j = 1,2, ...,p, that is 

'Gn(s) Gv 2 (s) . 

G'2l(s) 

. G, 


G(s) = 




Gpi{s) 


(«) 


Gpr( s) \ 


(1.18) 


Recall from basic circuits courses that while finding Gji(s) = Yj(s)/Ui(s ) all 
other system inputs except for Ui(s) must be set to zero. Similarly, the discrete¬ 
time transfer function of a multi-input multi-output, time invariant, system is 
given by 


G (~) = 


G\\(z) G\2 (z) 

. . . Glr(^) 

G 21 (^) 


Gji(z) 


G p i(z) 

• • • m 


pXr 


(1.19) 


More will be said about the transfer function of multi-input multi-output systems 
in Chapter 2. 


Since each entry in the matrix transfer functions given in (1.18) and (1.19) is a 
ratio of two polynomials with complex numbers, it is obvious that for higher-order 
systems the required calculations in the frequency domain are mathematically 
very involved so that the state space system representation is simpler than the 
corresponding frequency domain representation. In addition, since the state space 
method is based on linear algebra, for which numerous efficient mathematical 
methods have already been developed, the state space method is also more 
convenient from the computational point of view than the frequency domain 
method. However, the importance of the frequency domain representation lies in 
the simplicity of presenting some basic concepts, and hence it very often gives 
a better understanding of the actual physical phenomena occurring within the 
system. 
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1.4 Mathematical Modeling of Real Physical Systems 

Mathematical modeling of real-world physical systems is based on the application 
of known physical laws to the given systems, which leads to mathematical 
equations describing the behavior of systems under consideration. The equations 
obtained are either algebraic, ordinary differential or partial differential. Systems 
described by algebraic equations are of no interest for this course since they 
represent static phenomena. Dynamic systems mathematically described by 
partial differential equations are known as systems with distributed parameters. 
The study of distributed parameter systems is beyond the scope of this book. 
Thus, we will consider only systems described by ordinary differential equations. 
These systems are also known as systems with lumped (concentrated) parameters. 
Even the lumped parameter systems are, in general, too difficult from the point of 
view of solving the general control problem described in Section 1.1, so that we 
have to limit our attention to lumped parameter systems described by linear time 
invariant differential equations. Fortunately, many control systems do have this 
form. Even more, control systems described by nonlinear differential equations 
can very often be linearized in the neighborhood of their nominal (operating) 
trajectories and controls assuming that these quantities are known, which is very 
often the case, so that nonlinear systems can be studied as linear ones. The 
linearization procedure will be independently considered in Section 1.6. 

In the following the modeling procedure is demonstrated on a simple RLC 
electrical network given in Figure 1.4. Assume that the initial values for the 


h L R i 



Figure 1.4: An RLC network 
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inductor current and capacitor voltage are zero. Applying the basic circuit laws 
(Thomas and Rosa, 1994) for voltages and currents, we get 

e i[t) = L — + Ri'ii(t) + e o(f) (1.20) 


eo(/) = i? 2 * 2 (f) = ^7 j *3(r)dr =>• i 3 = 

o 

*l(f ) = *2(f ) + *3 (t) 

Using (1.21) in (1.22) produces 


hit) 


1 

R^2 


e o{t) + C 


de 0 {t) 

dt 


( 1 . 21 ) 

( 1 . 22 ) 


(1.23) 


Taking the derivative of (1.23) and combining (1.20) and (1.23), we obtain the 
desired second-order differential equation, which relates the input and output of 
the system, and represents a mathematical model of the circuit given in Figure 1.4 


d 2 eo(t) fL + R\R 2 C\ deo(t) f Ri R- 2 \ . . 

dt' 2 + V 1(21/’ ) dt + V R 2 LC ) e ° ( ) 


1 

LC 


e;{t) 


(1.24) 


Note that in this mathematical model r; (/) represents the system input and e 0 (t) 
is the system output. However, any of the currents and any of the voltages can 
play the roles of either input or output variables. 

Introducing the following change of variables 


a-T (t) 


A> 


dxi(t) 

dt 


d r o ( t ) 
dt 


z 2 (t) 


de 0 {t) 

dt 


z 2 (t) 


(1.25) 


U(t) = €i(t) 

y(t) = e 0 (t.) => y(t) = .ri (/) 

and combining it with (1.24) we get 


dx 2 {t) fL + RiR 2 C\ (Ri + R 2 \ t ^ 1 t ^ 

— + ( R 2 LC r m+ \~RtLC~ j' Tl * = LC nm <u6) 
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The first equation in (1.25) and equation (1.26) can be put into matrix form as 


X\ 


&2 m 



0 

.Ri -\-R 2 

r 2 lc 


L+R 1 R-,C 

R.LC 


,l’l 

X-2 


1 

.LC 


(1.27) 


The last eauation from (1.251. 


y(t) = [i 



(1.28) 


Equations (1.27) and (1.28) represent the state space form for the system whose 
mathematical model is given by (1.24). The corresponding state space matrices 
for this system are given by 


A = 


0 

1 

, B = 

' 0 ' 

^1 +^2 

L+R^R.C 

1 

L R.LC 

R.LC \ 


ILC J 


C = [1 0], D = 0 (1.29) 


The state space form of a system is not unique. Using another change of variables, 
we can get, for the same system, another state space form, which is demonstrated 
in Problem 1.1. 

The transfer function of this single-input single-output system is easily 
obtained by taking the Laplace transform of (1.24), which leads to 


G(s) 




s 2 + 


1 

LC 


L-\-RiR 2 C c 
R 2 LC * 


Ri ~\~R2 

r 2 lc 


(1.30) 


Note that a systematic approach for getting the state space form from 
differential (difference) equations will be given in detail in Chapter 3. In this 
chapter we present only the simplest cases. These cases are dictated by system 
physical structures described by a set of first- and second-order differential 
(difference) equations, which can be put in a straightforward way into matrix 
form, which in fact represents the desired state space form. 

Another example, which demonstrates how to get a mathematical model for 
a real physical system, is taken from mechanical engineering. 
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Figure 1.5: A translational mechanical system 


A translational mechanical system is represented in Figure 1.5. The following 
two equations of motion for this system can be written by using the basic laws 
of dynamics (Greenwood, 1988) 

F ' = "“if + Bl (ilr ^ ir ) + hin - m) <L31) 


and 


„ d 2 y 2 dy-2 

F-2 = m 2 —- + B 2 -j- + k 2 y2 ~ 
dt z dt 


~ % ) - - y -2 ) (1-32) 


dt 


It can be seen that this system has two inputs, iq and l >, and two outputs, tj\ 
and y >. The rearranged form of equations (1.31) and (1.32) is given by 


mi 


d 2 yi 

dt 2 


B 1 d ^ + k 1 y 1 - Bi^ - kiy 2 = F x 


(1.33) 


and 

— -Bi—^— k\yi + in 2 + [B\ + -B 2 )—+ (ki + k 2 )y 2 = F 2 (1-34) 

From equations (1.33) and (1.34) the state space form can be obtained easily by 
choosing the following state space variables 

dyi dy 2 

V 1? ^2 i ^3 £/2? ^4 


Ui = iq, u 2 = F 2 


(1.35) 
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The state space form of this two-input two-output system is given by 


and 


X\ 

X 2 

%3 

= 

'01 0 O' 

ki Bi ki Bi 

mi mi mi mi 

0 0 0 1 


X\ 

X 2 

%3 

+ 

' 0 O' 

— 0 

mi 

0 0 

X 4 _ 


ki Bi ki~\~k 2 Bi~\~B 2 

L m 2 m 2 m 2 m 2 -1 


x 4_ 


o 




Vi" 


'1 

0 

0 

O' 

.72. 


_0 

0 

1 

0_ 


X 4 


Ml 

.“2 . 

(1.36) 


(1.37) 


It is interesting to find the transfer function for this multi-input multi-output 
system. Taking the Laplace transforms of (1.33) and (1.34), and assuming that all 
initial conditions are equal to zero, we get the scalar transfer functions from each 
input to each output. This is obtained by keeping the input under consideration 
different from zero and setting the other one to zero, that is 


Gn(s) = 

G\2{s) = 

G' 2 l{s) = 

Gr' 2 2 (-S ) = 


Yi{s) 

UiG); | 


Lr 2 (,) = 0 


> |(-i 

U 2 (s) / | 

Y 2 (s)\ 


u-As)=o 




m-2S 2 

+ 

(Bi 

+ 

B 2 )s 

+ 

(h 

+ 

k 2 ) 

(I4S 4 

+ 

a 3 s 3 

+ 

a 2 s 2 

+ 

d\S 

+ 

a 0 



Bn 

‘ + ki 

1 




(I4S 4 

T 

a 3 s 3 

+ 

CI2S 2 

T 

d\S 

T 

ao 



Bn 

1 + ki 

L 




(I4S 4 

T 

a 3 s 3 

+ 

a 2 s 2 

T 

d\S 

T 

a 0 


ni\s 2 

+ 

B\s 

+ 

ki 



(I4S 4 

+ 

a 3 s 3 

+ 

CI2S 2 

+ 

d\S 

T 

a 0 


(1.38) 


a 4 = miiri2 

a 3 = 5i(mi + m 2 ) + miB 2 
a 2 = BiB 2 + k 1 (m 1 + m 2 ) + k 2 m 1 
(ii = k\B2 + &2-Bi 
a 0 = kxk2 


where 
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so that the system transfer function is given by 


G(.s) 


G\\(s) G'l 2 (.s) 

G 2 i(s) G 22 (.s) 


(1.39) 


Sometimes due to the complexity of dynamic systems it is not possible to 
establish mathematical relations describing the dynamical behavior of the systems 
under consideration. In those cases one has to use experimentation in order to 
get data that can be used in establishing some mathematical relations caused 
(induced) by system dynamics. The experimental way of getting system models 
is the subject of the area of control systems known as system identification. More 
about system identification can be found in Chapter 10. A classic textbook on 
system identification is given in the list of references (Ljung, 1987). 

The reader particularly interested in mathematical modelling of real physical 
systems is referred to Wellstead (1979) and Kecman (1988). 


1.5 Models of Some Control Systems 


In this section we consider the modeling of two common physical control systems. 
In that direction modeling of an inverted pendulum and a complex electrical 
network is presented. Mathematical models of a DC motor will be derived in 
Section 2.2.1. DC motors are integral parts of several control system schemes. 
Mathematical models of many other interesting control systems can be found in 
Wellstead (1979), Kecman (1988), and Dorf (1992). 

Inverted Pendulum 


The inverted pendulum is a familial - dynamic system used very often in 
textbooks (Kwakernaak and Sivan, 1972; Kamen, 1990). Here, we follow 
derivations of Kecman (1988) for an idealized pendulum of length / whose 
mass, mi, is concentrated at its end. It is assumed that a cart of mass m 2 
is subjected to an external force R | which, as a control variable, has to keep the 
pendulum upright. Cart displacement is denoted by x and pendulum displacement 
is represented by angle 9 (see Figure 1.6). 

Using elementary knowledge from college physics (see for example Serway 
1992), the equation of motion for translation in the direction of x axis is obtained 
by applying Newton’s law 


m 2 


d 2 x(t ) d 2 (x + / sin 9 ) 


dt 2 


mi 


dt 2 


= F 


(1.40) 
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Figure 1.6: Inverted pendulum 
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Equations (1.43) represent the desired mathematical model of an inverted pen¬ 
dulum. 

Complex Electrical Network 

Complex electrical networks are obtained by connecting basic electrical 
elements: resistors, inductors, capacitors (passive elements) and voltage and 
current sources (active elements). Of course, other electrical engineering elements 
like diodes and transistors can be present, but since in this course we study only 
the lineal - time invariant networks and since this textbook is intended for all 
engineering students we will limit our attention to basic electrical elements. The 
complexity of the network will be manifested by a large number of passive and 
active elements and large number of loops. Such a network is represented in 
Figure 1.7. 


Ri L\ L 2 L 2 R 4 



Figure 1.7: Complex electrical network 

In electrical networks composed of inductors and capacitors the total number 
of inductors and capacitors indicates the order of the dynamic system. In this 
particular' example, we have three inductors and two capacitors, i.e. the order 
of this dynamical system is n = 5. Having a dynamic system of order n = 5 
indicates that the required number of first-order differential equations to be set 
is five. If one sets more than n differential equations, for a system of order n, 
some of them are redundant. Redundant equations have to be eliminated since 
they do not carry any new information; they are just linear combinations of the 
remaining equations. Using basic laws for currents and voltages, we can easily 
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set up five first-order differential equations; for voltages around loops 

dii 

e; - L 1 — - Rill - v c = 0 
at 

~ + F(*2 - *3) = 0 


(1.44) 


and for currents 



+ i? 3 ( *3 - *4) - i? 2(*2 


* 3 ) = 0 


r , dv Cl 

*1 - *2 - Ci —— — 0 


dt 


i 4 ~ ( 2 


dv r 


dt 


= 0 


(1.45) 


We have set up five equations for four currents and two voltages. Current i \ can 
be eliminated by using the following algebraic voltage balance equation, which 
is valid for the last loop 


#3(43 - *4 ) = Rik + v C2 


(1.46) 


from which the cuiTent i \ can be expressed as 

R 3 . 1 

*4 = -*3 — - V, 

f ?3 + f ?4 f ?3 + f ?4 


(1.47) 


Replacing current *4 in (1.44) and (1.45) by the expression obtained in (1.47), 
the following five first-order differential equations are obtained 


dii 

dt 

di-2 

dt 


Ri . 1 1 

7" “ IT*'''+ 


R 2 ■ , R '2 ■ 1 

-~r i 2 + ^-* 3 + -j-v Cl 

T 2 T3 L 2 


d [3 

dt 


i ?2 • R2R3 + R2R4 + R3R4 ■ R3 

Tz 2 X 3 (jf? 3 + i? 4 ) * 3 “ £3(i23 + i24)'' 

dr’d 1 . 1 . 

FT = F * 1 ~ C~i 2 


dv r 


Rs 


*3 - 


d/ C' 2 (i (3 + i (4 ) C' 2 (i (3 + # 4 .) 


''c 2 


(1.48) 
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The matrix form of this system of first-order differential equations represents the 
system state space form. Take aq = ii,x - 2 = i-2,%3 = *3,aq = v Cl ,x$ = v C2 , 
u = <then the state space form is given by 


".i’r' 


X -2 


X3 

= 

X 4 


5 . 



Ri 

Ti 

0 

0 

X 

c\ 

0 


0 

_ Rj_ 

l 2 
Rl 

Rz 

_ x 

c\ 

0 


Ri. 

l 2 

« 33 

0 

R* 

T3 (Rz ~\~Rz ) 


_ x 

L 1 

x 

l 2 

0 

0 

0 


0 

0 

Rz 

L?j {Rz ~\~Rz ) 

0 

C' 2 {Rz-\-Rz) -I 



",1’1' 


r 1 -1 

L 1 


X 2 


0 



+ 

0 


X 4 


0 


. x 5. 


. 0 . 


(1.49) 


"'hr it n - — R 2 Rz-\-R 2 Ra-\-RzRa 

wheie a 33 - l^r.+r,) • 

Note that while modeling electrical networks, it is advisable to use the 
currents through inductors and the voltages on capacitors for the state space 
variables, which, in fact, is done in this example. 


We would like to point out that in this chapter we have presented only 
modeling and mathematical models for continuous-time real dynamic physical 
systems. The reason for this is twofold: (1) there are no real-world discrete¬ 
time physical dynamic systems; ( 2 ) discrete-time models obtained from social, 
economic, hydrological, and meteorological sciences are usually of no interest to 
control engineers. However, discrete-time models obtained by discretization of 
continuous-time systems will be treated in detail in Chapter 3. 


1.6 Linearization of Nonlinear Systems 

We have mentioned before that in this book we study only time invariant linear 
control systems and that the study of nonlinear control systems is rather difficult. 
However, in some cases it is possible to linearize nonlinear control systems 
and study them as linear ones. In this section we show how to perform 
linearization of control systems described by nonlinear differential equations. 
The procedure introduced is based on the Taylor series expansion and on the 
knowledge of nominal (operating) system trajectories and nominal control inputs. 
Readers particularly interested in the study of nonlinear systems are referred to 
a comprehensive book by Khalil (1992). 
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We will start with a simple scalar first-order nonlinear dynamic system 
represented by 


x{t) = f(x(t),u(t)) (1.50) 

Assume that under usual working circumstances this system operates along the 
trajectory x n (t) while it is driven by the control input u n (t). We call x n (t) and 
u n (t I* respectively, the nominal system trajectory and the nominal control input. 
On the nominal trajectory the following differential equation is satisfied 

.'iA(f) = f{x n {t),u n {t)) (1.51) 

Now assume that the motion of the nonlinear system (1.50) is in the neighborhood 
of the nominal system trajectory and that the distance from the nominal trajectory 
is small, that is 


x(t) = x n (t) + A x(t) (1-52) 

where A x(t) represents a small quantity. It is natural to assume that the system 
motion in close proximity to the nominal trajectory will be sustained by a control 
input which is obtained by adding a small quantity to the nominal control input, 
that is 


u(t) = u n (t) + Au(t) (1.53) 

For the system motion in close proximity to the nominal trajectory, from equations 
(1.50), (1.52), and (1.53), we have 

x n (l) + A x(t) = f(x n (t ) + A x(t),u n (t) + A u(t)) (1.54) 

Since A x(t) and A u(t) are small quantities, the right-hand side of (1.54) can 
be expanded into a Taylor series about the nominal trajectory and control, which 
produces 

Of Of 

X n{lj T Ax(t) — f (,T n {t ), U n (t)) ^ ( -I' n . u n )Ax(t) ^ ^ n . ll n )All(t) 

+ higher—order terms 


(1.55) 
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Using (1.51) and canceling higher-order terms (which contain very small quan¬ 
tities Ax 2 ,Au 2 ,AxAu,Ax 3 ,...), the following linear differential equation is 
obtained 

Ax(t) = -f-(x n , u n )Ax(t) + -f-(x n , u n )Au(t) (1.56) 


Ox 


du 


whose solution represents a valid approximation for A x(t). Note that the partial 
derivatives in the linearization procedure are evaluated at the nominal points. 
Introducing the notation 


a 0 {t) 


dj_ 

dx 






h _ df { v 


the lineal - system (1.56) can be represented as 


(1.57) 


A x(t) + a 0 (t)Ax(t) = b 0 (t)Au(t) (1.58) 


In general, the obtained linear system is time varying. Since in this course we 
study only time invariant systems, we will consider only those examples for 
which the linearization procedure produces time invariant systems. It remains to 
find the initial condition for the linearized system, which can be obtained from 
(1.52) as 


Ax(to) = x(t 0 ) - x n (t 0 ) (1.59) 

Similarly, we can linearize the second-order nonlinear dynamic system 

X = /( x,x,u,u) 

by assuming that 

x(t) = x n (t) + A x(t), x(t) = x n (t) + A x(t) 
u(t) = u n (t) + Au(/)| ii(t) = ii n (t) + A ii(t) 

and expanding 

x n + Ax = f( x n + Ax n ,x n + Ax,u n + Am, u n + Am) 

into a Taylor series about nominal points x n ,x n , u n , u n , which leads to 

Ax(t) + aiA.i’(/) + a 0 Ax(t) = biAii(t) + b 0 Au(t) (1.63) 


(1.60) 


(1.61) 


(1.62) 
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where the corresponding coefficients are evaluated at the nominal points as 


df . 

^1 — \%n i i i ^n) i 

a 0 = 

df , . 

i _ d f i i 

b 0 = 

d f i •, .. v 


The initial conditions for the second-order linearized system are easily obtained 
from (1.61) 


Ax(to) — - 1 (6)) % n{t o)i Ax(to) — ( 6)) &n{to) (1.65) 


Example 1.1: The mathematical model of a stick-balancing problem is given 
in Sontag (1990) by 

0 = sin 6 — u cos 9 = / (0, u} 


where u is the horizontal force of a finger and 9 represents the stick’s angular 
displacement from the vertical. This second-order dynamic system is linearized 
at the nominal points ( 9 n (t ) = 0 n (t) = 0 ,u n (t) = 0) by using formulas (1.64), 
which produces 


df n 

cii =-- = 0 , 

89 


a 0 = - 


dj_ 

89 


-(cos 9 -f u sin 9 )| 9„(o=o 


= -1 


(o=o 


/ d f n 

6 ' = m = °’ 


7 _ , df , 

On — I tt — I — 
8u 


— (cos 9) | = — 1 


1 Bn(t) = 0 


The linearized equation is given by 


9(t) — 9(t) = —u(t) 


Note that A0(t) = 0(t),Au(t) = u(t) since 0 n {t) = 0,u n (t) = 0. It is important 
to point out that the same linearized model could have been obtained by setting 
sin#(f) ~ 9(t ), cos 9(t) ~ 1, which is valid for small values of 0(t). 


o 

Of course, we can extend the presented linearization procedure to an //-order 
nonlinear dynamic system with one input and one output in a straightforward 
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way. However, for multi-input multi-output systems this procedure becomes 
cumbersome. Using the state space model, the linearization procedure for the 
multi-input multi-output case is quite simple. 

Consider now the general nonlinear dynamic control system in matrix form 
represented by 


^X(f ) = .'/'I XI / !. Ill / ! ! (1.66) 

where x(t), u(/), and T are, respectively, the n -dimensional state space vector, 
the r-dimensional control vector, and the n -dimensional vector function. Assume 
that the nominal (operating) system trajectory x n (/) is known and that the nominal 
control that keeps the system on the nominal trajectory is given by u n (t). Using 
the same logic as for the scalar case, we can assume that the actual system 
dynamics in the immediate proximity of the system nominal trajectories can be 
approximated by the first terms of the Taylor series. That is, starting with 

x(t) = x n (t) + A x(t), u(t) = u n (t) + Au(t) (1.67) 


and 

^x n (/) = 7(xj l (f), u^ff)) (1.68) 

we expand equation (1.66) as follows 


d d 

—x n + —Ax = /(x„ + Ax, u n + Au) 


= T(x„,u„) + 



Ax + 



Au f higher—order terms 

to 

in 

(1.69) 


Higher-order terms contain at least quadratic quantities of Ax and Au. Since 
Ax and Au are small their squares are even smaller, and hence the higher- 
order terms can be neglected. Using (1.67) and neglecting higher-order terms, 
an approximation is obtained 



Ax(f) + 

to 

in 



Au(f) 

in 

in 


(1.70) 
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Partial derivatives in (1.70) represent the Jacobian matrices given by 


d_r_ 

<9x 




i-TTl 

TTl 


&T, -i 



9a; i 

dF 2 

rJx 2 


* * * dx n 

a 


_ j^nXn _ 

dxi 


dT t 

dx n 

x„(0 




dxj 


■ .['1 


dZn. 

dT n 


dZn. 



■ dxi 

dx 2 


dx n ~ 


(1.71a) 


d_T_ 

<9u 


_ graXr _ 


'-nt') 


r 

TTl 


TTl 


c>u 2 


* * * du r 




a 




du r 



duj 


dF n 

'iZn. 


dT n , 


fJu 2 


du r 


|x„(() 


(1.71b) 


Note that the Jacobian matrices have to be evaluated at the nominal points, i.e. 
at x„(t) and u n (t). With this notation, the linearized system (1.70) has the form 

^Ax(t) = AAx(t) + BAu(t), Ax(/ 0 ) = x(/ 0 ) - x n (/ 0 ) (1-72) 

The output of a nonlinear system, in general, satisfies a nonlinear algebraic 
equation, that is 


y(t) = £(x(/),u(/)) 


(1.73) 


This equation can be also lineaiized by expanding its right-hand side into a Taylor 
series about nominal points x n (/) and u„(/). This leads to 


y n + Ay = t/(x n ,u n ) + 



Ax + 


+ higher—order terms 



Au 

to 

to 


(1.74) 


Note that y n cancels term (7(x„. u„ ). By neglecting higher-order terms in (1.74), 
the lineaiized part of the output equation is given by 


Ay (/) = CAx(<) + DAu(t) 


(1.75) 
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where the Jacobian matrices C and D satisfy 


(d_Q _\ 

V^X / ,*„(,) 


j-)pXr 



rdQ, 



ag i 

dxi 

dx 2 



ag 2 



dxi 


dg t 

* * * 



dxj 


dbl 




- dxi 

dx 2 


* * * dx n 


rdg 1 

9g, 


sg i 


c>u 2 


* * * du r 







dg t 

9it r 



duj 



dbl 



■ 

du 2 


9ti r 


(O 

(O 


(O 

to 


(1.76a) 


(1.76b) 


Example 1.2: Let a nonlinear system be represented by 

dx i 

—— = .l’i Sill X -2 + X- 2 U 

dt 


dx 2 
dt 


ie 


y = 2 x i x 2 + . 1*2 

Assume that the values for the system nominal trajectories and control are known 
and given by x lni .r 2 „. and u n . The linearized state space equation of the above 
nonlinear system is obtained as 


A.i’i(/) 


si n r 2) i 

•lira COS l 2n + «ra 

A.ii(f) 


%2n 

Ax 2 (t) _ 


_ e -1 ’ 2 " 

-hraf" 12 " 

Ax 2 (t) m 

1 



A y{t) = [2x 2n 


Arm T 2 .i*2,i 


Axi(t) 
Ax 2 (t) 


OA u(t) 


Having obtained the solution of this lineaiized system under the given control 
input A n{ I ), the corresponding approximation of the nonlinear system trajectories 

Axi(t) 

Ax 2 (/)_ 


x n (t.) + A x(f) = 


•i’i n( f) 

•r 2 n(f) 


o 
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Example 1.3: Consider the mathematical model of a single-link robotic 
manipulator with a flexible joint (Spong and Vidyasagar, 1989) 

I6\ + mgl sin 9\ + k{ 9\ — 0 3 ) = 0 

JO -2 - k{9 1 - 0 2 ) = u 

where 0 \ .0> are angular - positions, /,./ ar - e moments of inertia, m and / are, 
respectively, the link’s mass and length, and k is the link’s spring constant. 
Introducing the change of variables as 

■I’l = 01 , X 2 = 01 , X 3 = 0 2 , F4 = #2 

the manipulator’s state space nonlinear model equivalent to (1.66) is given by 

iq = X2 
mgl k 

X 2 = - —SU1X 1 - j(x 1 - x 3 ) 

X3 = X4 

k 1 

' 1>4 = 7 Xl ~ X3 + J u 

Take the nominal points as {xi n , X 2 n , x 3n ,X 4 n ,u n ), then the matrices A and B 
defined in (1.71) are given by 


0 

1 

0 

O' 


'O' 

k-\-mgl cos Xi n 

0 

k 

0 


0 

1 

I 

, B = 

0 

0 

0 

1 

0 

k 

L j 

0 

k 

J 

0 . 


r 

L j J 


In Spong (1995) the following numerical values are used for system parameters: 
mgl = 5, I = •/ = 1, k = 0.08. 

Assuming that the output variable is equal to the link’s angular - position, 
that is 


y = X! 


the matrices C and D, defined in (1.76), are given by 

C = [1 0 0 0], D = 0 


o 
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In the next example, we give state space matrices for two linearized models 
of an F-15 aircraft obtained by linearizing nonlinear equations for two sets of 
operating points. 

Example 1.4: F-15 Aircraft 

The longitudinal dynamics of an F-15 aircraft can be represented by a 
fourth-order mathematical model. For two operating conditions (subsonic and 
supersonic) two linear mathematical models have been derived (Brumbaugh, 
1994; Schomig et ah, 1995). The corresponding state space models are given by 


"i’i 1 [—0.00819 -25.70839 

x-2 _ -0.00019 -1.27626 

.7 3 “ 0.00069 1.02176 

_ 0 0 


0 -32.170951 pi-r 

1.0000 0 x -2 

— 2.40523 0 .x-3 

1.0000 0 . _.T4_ 


' -6.80939 
, -0.14968 

+ -14.06111 
0 

for subsonic flight conditions, and 


,i’i 
X 2 
X3 
X4 


71 -0.01172 -95.91071 0 -32.11294 

7 2 _ -0.00011 -1.87942 1.0000 0 

7 3 “ 0.00056 -3.61627 -3.44478 0 

7 4 0 0 1.0000 0 



-25.40405 
-0.22042 
-53.42460 
0 

for supersonic flight conditions. 

Model derivations are beyond the scope of this book. The state space 
variables represent: x\{t )—velocity in feet per second, x 2 {t )—angle of attack 
in radians, ./■;;(/)—pitch rate in radians per second, and ./■ ](/)—pitch attitude in 
radians. The control input u(t) represents the elevator control in radians. 




o 
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Linearization of an Inverted Pendulum 

Sometimes is not necessary to go through the entire linearization procedure. 
It is possible to simplify and linearize mathematical equations describing a given 
dynamic system by using simple mathematics. This will be demonstrated on an 
example of the inverted pendulum considered in Section 1.5. A linearized model 
of the inverted pendulum can be obtained from equations (1.43) by assuming that 
in the normal operating position (pendulum in an upright position) the pendulum 
displacement 9 is very small so that the following approximations are valid 
sin 9(t) ss 9(t), cos 9(t) ss 1, 9(t)(d9(t)/dt) 2 ss 0, Vt Then, from (1.43), the 
linearized model of the inverted pendulum is obtained as 


d 2 x d 2 0 

[,n 1 + 1 n 2 )- m + ,n 1 l li¥ = F 


d 2 x d 2 0 
dt 2 + dt 2 


g® 


(1.77) 


This model can easily be put in the state space form equivalent to (1.14) by 
introducing the following change of variables 


Ti = x =>• ,i’i = x-2 
dx 

X-2 = — = X 

dt 

, 1' 3 = 6 =>- . i ’3 = , 1’4 

110 
= Tt 
u = F 

With this change of variables equations (1.77) imply 

mi g 1 

To =- x 3 -\ - u 

ni 2 ni 2 

(m 1 + m 2 )g 1 

T’4 = ----—T’3-7 U 


m 2 I 


m 2 l 


(1.78) 


(1.79) 


From (1.78) and (1.79) the state space form of the inverted pendulum is given by 


X\ 


'0 

1 

0 

O' 


X\ 


' 0 ' 

Xo 


0 

0 

_ rn l g 

0 


Xo 


i 


— 

m 2 



+ 

m 2 



0 

0 

0 

1 




0 

X4_ 


.0 

0 

(m 1 +m 2 )g 
m 2 l 

0. 


_ x 4 _ 


_r_ 

m 2 l - 


u 


(1.80) 
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The output equation can be chosen such that it gives information about the call’s 
horizontal displacement aq and the link’s angular position a 3 . In that case a 
possible choice for the output equation is 

y=[ 1 0 1 0 ]x (1.81) 


The same state variables will appeal - directly on the output if the following output 
equations are used 


y = 


l 

o 


0 0 0 
0 1 0 


(1.82) 


y = 


o o 
1 0 


1 0 
0 0 


X 


(1.83) 


The main difference between (1.81) and (1.82)—(1.83) is that in output equations 
(1.82) and (1.83) we have two channels, each of which independently produces 
information about the particular state variable. 

Finally, we would like to point out that the SIMULINK package is very 
convenient for simulation of nonlinear systems. It can also be used to obtain 
linearized models of nonlinear systems around given operating points (nominal 
system trajectories and controls). 


1.7 MATLAB Computer Analysis and Design 

MATLAB is a very advanced and reliable computer package, which can be used 
for computer-aided control system analysis and design. In addition to handling 
standard linear algebra problems, it has several specialized control theory and 
application toolboxes. One of them, the CONTROL toolbox, will be extensively 
used in this book. At some places in the book we also refer to the SIMULINK 
package, which is very convenient for simulation (finding system responses due to 
given inputs) of linear and nonlinear systems. MATLAB is user friendly. It takes 
only a few hours to master all of its functions. MATLAB will help students obtain 
a deeper understanding of the main control theory concepts and techniques, and 
to study higher-order real physical control systems, which would be impossible 
using only pencil and paper. Many MATLAB problems, laboratory experiments, 
and case studies will be encountered in this book. More about MATLAB and its 
CONTROL toolbox can be found in Appendix D. 
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1.8 Book Organization 

This book has two parts: Part I, titled Methods and Concepts , Chapters 2-5, 
and Part II, titled Analysis and Design , Chapters 6-9. In the first part of the 
book, in Chapters 2 and 3, two main techniques of control theory—the transfer 
function approach and the state space method—are discussed in detail. In Chapter 
4 we consider the concepts of system controllability and observability, and in 
Chapter 5 the stability concept of time-invariant continuous and discrete systems 
is presented. 

In this introductory chapter we have defined the general control problem. 
The main control system characteristics and control objectives will be presented 
in the second part of this book stalling with Chapter 6. In Chapter 6, the control 
system specifications relating to a system’s transient behavior and steady state 
properties will be considered. Since the emphasis in this book is on the time 
domain controller design (based on the root locus technique), the corresponding 
technique is presented in detail in Chapter 7. Design of controllers that solve 
specific control problems will be presented in Chapters 8 and 9. In Chapter 10 we 
give an overview of modern control theory, which can serve as an introduction 
for further studies of control theory and its applications. The presentation of 
Chapter 10 to undergraduate students can either be completely omitted or even 
expanded at schools that have strong programs in control systems. 
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1.10 Problems 

1.1 Find the state space form for the electrical circuit whose mathematical model 
is given in (1.20)—(1.22) by taking for the state space variables the input 
current and output voltage, i.e. by choosing ,i’i = i\,x - 2 = e 0 . In addition, 
take y = x - 2 and u = e% 



INTRODUCTION 


35 


1.2 Find a mathematical model for the inverted pendulum, assuming that its 
mass is concentrated at the center of gravity, and linearize the obtained 
nonlinear system at 0 n (t) = 0 n (t) = x n (t) = x n (t) = 0 (Kwakernaak and 
Sivan, 1972; Lewis, 1992). Compare the linear model obtained with model 
(1.80) derived under the assumption that the pendulum mass is concentrated 
at its end. 

1.3 Verify the expressions given in (1.38) for the transfer function of a two-input 
two-output translational mechanical system. 

1.4 Find the mathematical model, transfer function, and state space form of the 
electrical network given in Figure 1.8. 



Figure 1.8: An electrical network 

1.5 Linearize a scalar system represented by the first-order differential equation 

= x (t)u(t)e~ u ^\ ,r(0) = 0.9 

at a nominal point given by (x n (t),u n (t)) = (1,0). 

1.6 Consider a nonlinear continuous-time system given by 

d 2 x(t) dx(t) dx(0) 

~ = - 2 ^“ cos + u(t))x(t)+ 1, .r(0) = 1.1, = 0.1 

Derive its linearized equation with respect to a nominal point defined 
by (x n (t),u n (t)) = (1,0). Find the linearized system response due to 
A u(t) = e~ 2t . 
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1.7 For a nonlinear system 

d^(t)_ + + (x _ u(t )) x (t) = u 2 (t) + 1, .r(0) = 0, = 1 

find the nominal system response on the nominal system trajectory defined 
by u n (t) = 1, subject to ,xq(0) = 0 and dx n (0)/dt = 1.1. Find the 
linearized state space equation and its initial conditions. 

1.8 The mathematical model of a simple pendulum is given by (see for example 
Kamen, 1990) 

d'20 

I— + mgl sin 6 = lu(t), 9(t 0 ) = 9 0 , 9(t 0 ) = 

where I is the moment of inertia, /, m are pendulum length and mass, re¬ 
spectively, and u(t) is an external tangential force. Assume that 9 n (t) = 0, 
u n (t) = 0,9 n (t o ) = 0 ,9 n (t 0 ) = 0, and are small. Find the linearized 
equation for this pendulum by using formulas (1.64). Determine the initial 
conditions. 

1.9 Linearize the given system at a nominal point (i „ • ) = (0,1,1) 


iq = xix -2 — sin.i’i 
x -2 = 1 — 3x 2 e~ Xl 

■i '3 = X1X2X3 

1.10 Linearize a nonlinear control system represented by 

iq = rrln.xq + .xqe -11 
x 2 = .i’i sin u — sin .xq 
y = sin.xq 

Assume that xi n ,X 2 n , and u n are known. Find the transfer function of the 
linearized model obtained. 

1.11 Linearize the Volterra predator-prey mathematical model 


iq = — ,xq + .xq.i ’2 
i'2 = X 2 ~ X\X2 

at a nominal point given by (./■ 1 „. ./q,, ) = (0,0). 
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1.12 A linearized model of a single-link manipulator with a flexible joint is given 
by (Spong and Vidyasagar, 1989) 

JiOi + BA + m - 0 m ) = 0 

~\~ B m 6 m k( 6i 0m) — ) 

where Ji,-J m are moments of inertia, Bi,B m are damping factors, k is 
the spring constant, u(t) is the input torque, and 0 m (t.),0[(t) are angular 
positions. Write the state space form for this manipulator by taking the 
following change of variables: .iq = d/. ,r> = 0i,x 3 = 0 m ,X 4 = 0 m . 

Remark: Note that the SIMULINK package can be used for linearization of 
nonlinear systems. Students may check most of the linearization problems given 
in this section by using the linmod function of SIMULINK. 



